Converting a series in A to a series in A 
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We introduce a transformation for converting a series in a parameter, A, to a series in the inverse of 
the parameter A -1 . By applying the transform on simple examples, it becomes apparent that there 
exist relations between convergent and divergent series, and also between large- and small-coupling 
expansions. The method is also applied to the divergent series expansion of Euler-Heisenberg- 
Schwinger result for the one-loop effective action for constant background magnetic (or electric) 
field. The transform may help us gain some insight about the nature of both divergent (Borel or 
non-Borel summable series) and convergent series and their relationship, and how both could be 
used for analytical and numerical calculations. 
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PACS numbers: 11.15Bt, ll.lOJj, 11.15Me, 0.23.Lt 

INTRODUCTION 

Much of our understanding of physics and mathematics 
relies on our ability not only to find relationships between 
various quantities through equations (such as differential 
equations), but also their solutions. In many cases one 
resorts to assuming that the solution can be expressed as 
a series expansion ^c(n)A™ in some parameter A. De- 
pending on the nature of the problem at hand, a first 
attempt at a series expansion might yield a convergent se- 
ries, whereas in other situations, a divergent series. One 
finds in many instances that a divergent series arises, e.g. 
in perturbative expansions in field theories QLQQ 

Convergent series are generally well behaved. Diver- 
gent series, on the other hand, require more careful atten- 
tion to elucidate their meaning. In many scenarios, the 
information one obtains, such as an expression for the co- 
efficient function c(n), is very useful and various schemes 
has been proposed that allow to one assign values to di- 
vergent series. Many analytical and numerical methods 
have been put forward, e.g. Borel summation, Pade and 
Borel-Pade approximants, large order perturbation the- 
ory (see and references therein) , variational perturba- 
tion theory (|5| and refs. therein), sequence transforma- 
tions (|(| and refs. therein). 

Another, but little known, numerical procedure is that 
of Mellin-Barnes regularisation originally proposed by 
Kowalenko et al 0, H > which can be used to obtain 
finite numerical values of a divergent (or convergent) se- 
ries, for certain values of complex A. This method is an 
exact procedure in the sense that no approximations are 
made, and it yields the same numerical values as those 
obtained by Borel summation for Borel summable series. 

Further investigation of the numerical Mellin-Barnes 
regularisation method, particularly with coefficient func- 
tions c(n) containing T-functions, and extending A to the 
entire complex plane (i.e. including non-Borel summable 



series), revealed an analytical relationship between series 
expansions in A to those in A -1 . This is in contrast to 
the basic Borel transformation, in which case, before and 
after the transform one still has a series in A. 

With this new insight we show how one can trans- 
form a series in ^c(n)A™ to a series in ^ d(n)\~ n . It 
then becomes apparent that the transform makes no dis- 
tinction between Borel (alternating) and non-Borel (non- 
alternating) summable series. In fact, in some cases it 
converts one to the other, and also converts a divergent 
series to a convergent series. As we shall see, depending 
on the form of the coefficient function c(n) , a single series 
in A can yield more than one series in A -1 . 

Many functions have both convergent and divergent 
series, and integral, representations. Our ability to un- 
derstand the solution to a problem is only as good as our 
ability to extract information about the coefficient func- 
tion c(n). An advantage of the proposed transformation 
is that given c(n), we can look at both the series in A 
and in A -1 . If one has the exact expression for c(n), be 
it as a convergent or divergent series, then it can be re- 
garded as the solution to the problem. If we only have 
an approximate expression for c(n), it is hoped that con- 
verting ^c(n)A™ to ^G?(n)A~™ might allow us to gain 
more understanding of the solution to the problem. 

THE TRANSFORM 

Consider a sum of a finite number of terms 

N 



S{x, N) = J2 c ( n K' 



(1) 



n=0 



'Electronic address: a.rawlinson@physics.unimelb.edu.au 



for any (meromorphic) coefficient function c(n). Using 
Cauchy's integral theorem, we can express this as 

S(x,N)=I dt f^' EE / (2) 
Jc(N) exp(2TTit) - 1 J C{N) 

where we have introduced a short hand notation for the 
integral. It is assumed that the only poles within the 



2 



contour C(N) are the simple poles (of residue l/(27ri)) 
of l/(exp(27r«i) - 1) at t = 0, 1, 2, . . . , N. (Of course, if 
there are poles due to c(t), their residues are included). 

If we express the closed contour C(N) = L(N) + L 
where L intersects the Re t axis somewhere at — 1 < a < 
0, and L(N) intersects the Ret axis between N and N+l 
as shown in Fig.^ the nature of c(t)x*/ (exp(27rii) — 1) in 
the complex i-plane may allow one to adjust the angles 
< 0i < 7r and — 7r < 02 < of L so that 



for t 
and t 



lim - 
fi^oo cxp(27rii) — 1 

aexp(z7r) + i?exp(i0i) 

aexp(i7r) + i?exp(i0 2 ) 







(3) 



where we note the angles are centred at t = (— a, 0) = 
aexp(i7r), rather than at the origin f = (0,0). If these 
limits can be achieved, the integral along L is indepen- 
dent of N. This is possible if c(t) contains T-functions. 
See Fig. El for an example, where the coefficient function 
is given by (|27|l . Note that the choice of contours L and 
L(N) will depend on both x and c(t), since for a given 
c(t) it is possible one will need a different contour for 
x = 2, say, compared to x = — 2. 
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FIG. 1: The complex i-plane showing the regions enclosed by 
the contours C(N) = L{N) +L and C"(-M) = L'(-M) - L. 



Using the short hand notation in |J5J, we can express 
the loop integration in as 



TV 



C(N) JL(N) JL 

= 2m^Res{C{N)) = ^c{n)x n (4) 



where Res(C(N)) denotes the sum of the residues of 
the poles within C(N), where the last equality is based 
on the assumption that c(t) has no poles within C(N). 

Now consider another contour C"(— M) = L'(—M)—L, 
where L'(—M) intersects the Ret axis somewhere be- 
tween —M—l and —M, as shown in Fig. 2] and we note 



the sign in front of L is changed because of the counter- 
clockwise convention of Cauchy's intergral theorem. One 
finds that 



/ = / -/ =2wiy2Res(C'(-M)) 

Jc(-M) JL'(-M) JL ' 



'L'(-Af) 



(5) 



where Res(C'(—M)) is the sum of the residues of the 
poles within C'(—M). Thus, combining J3J and JSJ, we 
find that 



/' = 5>(n)*»- / 
Jl „ =0 Jl{ 



L(N) 

2iri^Res{C'{-M)) 



(6) 



L'(-M) 



This is the main result of the paper, and is applicable 
to cases where c(n) represents either a convergent or di- 
vergent series. One can use this result for analytical and 
numerical studies of a series. In numerical studies, if 
the contour L in Fig. is chosen with 0i = ir/2 and 
02 = — 7r/2, the first line of I© gives the Mellin-Barnes 
regularisation 

In general, one often finds that if J2n=o c i n ) xn is a 
divergent series, then Res(C'(—M)) is a convergent 
series, and the nature of the integrand in ijSJ) is such that 
as M — > 00 then j L t/ M \ — > 0, so 



N 

£< 

n=0 



(n)x n - [ = -2m Res(C (-00)) = f 
Jl(n) Jl 



(7) 

Therefore we find that — 2ni Res(C'{— 00)) is finite 
and can be regarded as the convergent series representa- 
tion of X)^=o c ( n ) a; "' an d a l so °f Jl- Observe that since 
Yln=o c { n ) x71 is a divergent series, it grows with TV, but so 

does Ji(iv)' in such a wa y that Sn=o C ( n ) xn -J L (N) = Jl 
is finite and independent of N, i.e. 



c(n)cc" 



E 

N' 

71=0 



L(N) 



L(N') 



(8) 



SOME EXAMPLES 

Let consider an example with a series having a coeffi- 
cient function of c(t) = T(t + a) where, without loss of 
generality, we take a > and a non-integer. Subsituting 
this into 10 we see that the only poles within C(N) are 
simple poles at t = 0, 1, 2, . . . due to l/(exp(27ri£) — 1), 
and their residues sum up to yield the divergent series 



TV 



2-xi ^2 Res(C{N)) = J2 T ( n + a ) xn - 



(9) 
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Within C'(—M) we find two sets of simple poles. One 
set at t = —1,-2,-3, . . . from l/(exp(27ri£) — 1) where 
the residue of the pole at t — — m is 



x m T(a — m) 



(— x) m sin(a7r)r(m + a — 1) 



(10) 



where we have used the T-function reflection formula, 
r(x)r(l — x) — 7r/ sin(7ra). Summing these residues, we 
see that it is a convergent series and therefore we can let 
M -> oo :- 



E 



— ' (— x) m sin(a7r)r(m + a — 1) 



7re 



-l/x 



(—x) a sin(a7r) 



1 



T(l-a,-l/x) 
r(i -a) 



(11) 



where T(y, z) is the incomplete T-function. This series, 
arising from the poles of l/(exp(27rii) — 1), can be re- 
garded as the counterpart of Jj5J where we take the sum 
over negative, instead of positive, integers. 

The other set of poles at t = —a, —a — 1, —a — 2, . . . are 
from r(t + a), and the residue of the pole at t — —a — m 
is 



(12) 



{-x) a+m sin(o!7r)r(m + 1) 



and summing up these residues gives a convergent series 
as well :- 



E 



-n (~x) a+m sin(o!7r)r(m + 1) {-x) a sin(a7r) 



(13) 



m— 



which will cancel the first term in the brackets in (|ll|l . 
This series arises from the poles in the coefficient function 
c(t) = r(i + a). Therefore, the total sum of residues 
within C'(— co) is 



2m^Res{C'{-oo)) = 



-e- l / x T{a)T{\ - a, -l/x) 



(14) 

where we have used the T-function reflection formula 
again. We thus see the relation of the divergent series © 
and the convergent series (|14|l via ©. 

If we now consider a to be an integer in JHJ, say, a = 1, 
we get the so-called 'paradigm divegerent series' :— 



N 

E 

ra=0 



r(n + l)x n . 



(15) 



Inserting this into JSJl, we have simple poles at t = 
0,1,2,... within C(N), and the sum of their residues 
yields the series (I15|) . Since T(t + 1) has poles at 
t = —1,-2,-3,..., we now have double poles within 
C'(—M), where the residue of the pole at t = —m is 



-(log(— x) + ip(m)) 
(-x) m T(m) 



(16) 



where ip(n) is the digamma function. Since the sum of 
these residues is a convergent series, we can take M — » oo, 
and obtain 



2m ^Res(C (-oo)) = 



-(log(-ar) +ip(m)) 



(~x) rn T(m) 



-VaT(0,-l/x) 



(17) 



Note that this is the same as (|14[l with a = 1. Notice the 
presence of log(— x) in each residue term in H16|l. which 
can pulled out of the summation in (T7J|. However, the 
summation of the i/j(m) terms yields a log(— l/x) term, 
such that it precisely cancels the log(— x), resulting in 
the log(— x)-free expression in Ijl7|l . This is an indication 
whereby if one finds logarithmic terms in a series expan- 
sion, that upon summation of the series such terms might 
vanish. 

Thus we see some of main features of © :- a series 
in x n is converted to a series in (—l/x) n , i.e. a non- 
alternating series is converted to an alternating series 
in the inverse of the expansion parameter. If double, 
or higher order, poles are encountered, derivatives (wrt 
t) of the coefficient function, c(t), and power term, x , 
yield derivatives of T-functions and log(x) terms respec- 
tively. We note that if x < 0, then l|17l) is real, whereas 
if x > 0, a non-zero imaginary part arises of magnitude 
— (tt/x) cxp(— l/x). 

If we consider the application of Borel summation on 
the 'paradigm divergent series' (PDS), (fH)|) . :- 



(18) 



one inserts 



p{x)pds = ^T(n+l)x", 



/•CO 

r(n + l)= / dse~ s x r - 
Jo 



(19) 



into Ijl8|l . interchanges the integral and summation signs, 
and then does the summation to yield 



p{x) 



Borel 



ds ■ 



1 — XS 



(20) 



If x < 0, we have an alternating series, and it is asserted 
that the 'finite value' of the divergent series (|18|l is 

P(x)pDs\ finite = p(x)Borel = ~Eq. (JT7J| X < (21) 

and yields the same results (with a minus sign) as in Ijl7|l . 
In this case the series is Borel summable. 

If, however, x > 0, we have a non-alternating series and 
there is now a pole on the line of integration at s = l/x 
in (|20|l . A potential ambiguity now arises as one would 
need a prescription (e.g. principal parts) to handle the 
pole. Such non-alternating series are regarded as being 



4 



'non-Borel summable'. These issues do not arise when 
one uses the transform since it can be used in ex- 
actly the same way to evaluate non-alternating series as 
for alternating series (with appropriate contours for each 
case). 

Let us now consider the convergent series representa- 
tion of the exponential function exp(a;) = J2^=o x n /n\. 
Putting c(t) = l/T(t + 1) in ©, then within C(N), the 
only poles are at t = 0,1,2,... due to 1 / (exp(2irit) — 
1), and the sum of the residues gives Xln=o a; "/ n ' 
and since this is a convergent series, we can take 
limAT— »oo Jl(n) 0- However, while 1/T(t+ 1) has ze- 
roes when t = —1, —2, —3, . . . (because L(z) has poles of 
residue (— l)"/n! at z = —n = 0, —1, —2, . . .), there are 
no poles within C'(—M) since the poles of 1 / (exp(2irit) — 
1) are cancelled by the zeroes of 1/F(i + 1). This means 
that exp(x) has the interesting feature that it has no di- 
vergent series representation. Therefore we find from © 

oo Tt r r 
27riVi*e S (C(oo))=exp(x) = V ^- = / = / 

n=0 n - J L JL'(-M) 

(22) 

where the last term is in fact independent of M. 

Generally it is the combination c(n)x n that determines 
whether the series is convergent or divergent, e.g. the ge- 
ometric series ^ x n has a finite radius of convergence - it 
converges if \x\ < 1, and diverges for x > 1. However, for 
many series it is primarily the nature of the coefficient 
function c(n), rather than the power term x n that de- 
termines their convergent or divergent nature. Knowing 
c(n) for some function, and if it gives rise to a divergent 
series, one can use this technique to obtain a convergent 
series representation of that function. It is also worth 
pointing out that given a convergent series representa- 
tion of a function, one can use this method to obtain its 
divergent series representation, if such exists. 

Let us look at an example of how a divergent series 
arises, and then try to gain an understanding of how 
such objects might be handled. Consider the following 
integral, 

r°° 2 4 e 1 /^) 
Z(g) = J o dxe-*s* = —=-Ki(l/(8g)) (23) 

where g > and K v {x) is the modified Bessel function. 
If we express exp(— x 2 ) as a power series, interchange the 
integration and summation, we find that 



n=Q vy 



T(n+1) 



(24) 



If, instead, we express exp(— gx ) as a power series and 
do the same operations, we see that 



"Z(9)" 



J> r(n+l) 

n=0 V ; 



(25) 



Using Stirling's formula F(n + 1) w exp(nlnn — n) 

for large n, we see that for large enough n, T(n) grows 
faster than x n for all x. If we take (12511 literally we see 
that the series diverges as we add more and more terms, 
which is why quotes are used in l(23|) . i.e. "Z(s)" ^ 
Z{g). Thus the 'small-g' expanson l(23|) gives a divergent 
series, while the 'large-g' expansion (|24l) is a convergent 
series. These are different series representations of the 
same function. 

The origin of the divergent behaviour in (|25(l is due to 
the fact that we have been 'careless' in interchanging the 
order of integration and summation (of an infinite num- 
ber of terms) in deriving Ij25(l . However, the coefficient 
function of the divergent series H25JI 



r(2r 



r(n + i) 



(26) 



does contain all the information about the function Z{g). 
To see this, consider the sum of a finite number of terms 

of (El 



n 



n=0 



r(n + l) 



N 

E 

n=0 



(-4g) n - 



T n 



Tin 



2V27rL(n+ 1) 

(27) 

where we have used L(2n) = 4"T(n)r(n + 1/2)/ '(2-v/tF) 
(E3 8.335). Inserting the coefficient function into © we 
find that the residues of the poles within C(N) gives 12711 . 
It is well known that for small g we can sum a few terms 
to give approximate (and very accurate) results for Z(g) 
- this is the essence of asymptotic expansions. But if we 
consider 



?i=0 



l)r(r 



(It 



2V27rL(n + 1) 



(28) 



L(N) 



2 v / 2rf\i + l)(exp(27ri<) - 1) 



it gives the same numerical values as the integral form of 
Z(g) for Re(g) > 0. Figure [5] gives a density plot of the 
real part of the integrand in l|29() for g = —2.5 showing 
the oscillatory nature, and also a suitable contour for L 
is given by the dark line. A density plot of the imaginary 
part of the integrand is similar. 

Within C'(—M) there are two sets of poles, each 
set giving a convergent series. Letting M — * 00, 
we find the residues of the set 'A' of poles at t = 
-1/4,-5/4,-9/4,... are 

Z A (J] y -V?(4g)-"r(n+l/4) 

l9) ~^ 4 ff V4 r( „ + l/2)F(n+l)' (29j 

while the set 'B' of poles of at t = -3/4,-7/4,-11/4, . . . 
gives residues 



Z' 



0F(4ff)~"F(n + 3/4) 
8 ff 3 / 4 F(n + l)F(n + 3/2) 



(30) 



5 




> o.i 



I 



< -0.1 
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FIG. 3: The imaginary parts of Z A (g), Z B (g) and Z A (g) + 
Z B {g) given by (|33J, for g < 0. 



FIG. 2: A density plot of the real part of ((-g)T(2t + 
l/2))/(2r(t + l)(exp(27rit) - 1)) for g = -2.5. Light and 
dark bands show the oscillatory nature of the function. The 
dark line is a suitable contour for L. 



and 

Z(g) = Z A (g) + Z B (g). (31) 
If g is real, the sums can be converted to Bessel functions, 



z A {g) 



z B (g) 



-ire 



V(8ff) 



4V¥ 

P 1/(8<J) 



4V¥ 



-A(l/(8ff)). 



(32) 



For the case of g < 0, we can express the imaginary part 
of Z(g) as 



Im{Z A - B (g)) 
Im{Z A (g)) 



(Z A > B (g + ie)-Z A > B (g-ie)) 



-l/(8|ff|) 



Jm(Z s ( 5 )) = 



-l/(8|ffl) 



-/_i(l/(8bl)) 



1.91 

7-W 



-Z 4 (l/(8| 5 |)) 



/m(Z( ff )) = Im(Z A (g))+Im(Z B (g)) 
Tre-i/Cslsl) 



-^i(V(8|5l)). 



(33) 



Figure shows a plot of the imaginary parts. 

ONE LOOP QED EFFECTIVE ACTION 

The Euler-Heisenberg-Schwinger one-loop QED effec- 
tive action [ToL HH IT^| , for a constant background mag- 
netic field B can be expressed as a proper time integral 



e 2 L 3 TB 2 
S = 5-5 — S B {gB) 



(34) 



where L 3 T is the space-time volume factor, and 



Sb {gi 



ds 



coth s 



exp 



<9b 



(35) 



where, <?s = e 2 B 2 /m i . If we expand coth(s) as a power 
series about s — we obtain the divergent series (and 
hence the use of quotes) :- 



"<Sb(.9b)" = 16 #b 

with coefficient function 
c(n) 



n=0 



c(«)5s 



(36) 



\y+H n \B 2n+i \ 



(2?i + 4)(2n + 3)(2n + 2) 

-C(-2n-3)(-4)" 
2(2n + 3)(n + l)cos(nTr) 



(37) 



and i?2n are the Bernoulli numbers, whose magnitudes 
are given by 



I -B 2 



(38) 



while £(s) is the Riemann (^-function. We have used the 
following relation between the T- and functions (0] 
9.535) :- 



T(z)((z) 



2cos(z?r/2) 



(39) 



Notice in (|37|l we have not written 1/ cos(n7r) = (—1)™ for 
n integer as we wish to maintain the analytical structure 
of (|39() for complex z. 

If one replaces l/cos(n7r) with (—1)" (even though 
they give the same divergent series), which will cancel the 
other (—1)™ term in [j3*7)l . the pole/singularity structure 
at t = —1, —3/2, —2 will give different residues from the 
simple (i.e. non-leading singularities) poles at l/(t+ 1), 
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l/(t + 3/2) and l/(t + 2). The resolution of ambiguities 
of this nature depends on the problem at hand. In this 
case, the term 1/ cos(z7r) is not the same as (— l) 2 for 
complex z. 

If we insert the coefficient function 1(3 7|) into 10 we 
find within the contour C(N) there are simple poles at 
t = 0, 1, 2, . . . whose residues give (@UJ>- Thus 



N 



Sb(9b) 



16 5 sv- C(-2r i -3)(-4 ffB ) J 



E 



^ (2n + 3)(n + l)cos(mr) y i(JV) 

(40) 

Within C'(—M) we have simple and double poles 
at t = —1,-3/2 and —2; two sets of simple poles, 
one set at t = —3, —4, —5, . . . and another set at t = 
—5/2,-7/2,-9/2,.... The residues of all poles within 
C'{—M) are related to Sb{9b) through 

S b {9b) = -2 m y2 Res(C'(-M))+ f . (41) 

J L'(-M) 

Letting M — » oo, we find the closed form 

Sb(9b) = i--^-4C'f-l,l + -^) 
3 HlB V 2 V9bJ 



\ + ^- + T7=) 1 °^9b) (42) 

6 4 5B 2^i y 



where C'( s j a ) is the derivative of the Hurwitz ^-function, 



C'(s,a) = — C(s,a). 
as 



(43) 



Numerical values for 1)42(1 can be obtained using Mathe- 
matica 0], where C'(~ L a) is expressed in terms of the 
Barnes' G-function |l5j. 

If instead one has constant electric fields, the one-loop 
effective action is 



e 2 L 3 TE 2 
S = — Se(9e) 



(44) 



where §e — e 2 E 2 /m 4 , and 

Se{9e) = -S b {-9e)- (45) 
One finds very good numerical agreement of 
Im(SE{gE)) = Im(—SB(—9E)) 

-wk\ 

(46) 



1 X - 1 

fe=i 



DISCUSSION AND CONCLUSION 

The above examples illustrate how various series ex- 
pansions of a function are related. If a problem is well 
posed, one would hope to obtain well defined results. If 



a divergent series in x arises, the starting point is to uti- 
lize the information contained in the coefficient function 
c(t). By considering a sum of a finite number of terms, 
i.e. a partial sum, and expressing it as a Cauchy inte- 
gral representation via one then attempts to find a 
contour C(N) that can be split into two parts L(N) and 
L as in Fig. If a suitable contour can be found, the 
integral along L(N) is subtracted from the partial sum, 
the result of which gives the integral along L. This in 
turn is related to the corresponding (convergent) series 
in x^ 1 obtained from the same coefficient function c(t). 

Equation ©, relating a series in i to a series in a;" 1 
is an exact relation. We also see there is no essential dif- 
ference between alternating (Borcl summable) and non- 
alternating (non-Borel summable) series, in contrast to 
the situation where a prescription is required to handle 
the potential ambiguity in the Borel summation proce- 
dure for non-alternating series (e.g. for the case of x > 
in (25)) )■ If one has a situation where there are several 
T- functions in the numerator of c(t) , then the Borel sum- 
mation procedure would result in each T-function being 
expressed in the integral form (|19[1 . leading to a multi- 
dimensional integral version of (|2t)|l . An advantage of 
using jnj is that it is a one-dimensional integral for arbi- 
trary c(t). 

If there are T- functions in the coefficient function in the 
form of r(n + a) where n is a positive integer, and a not 
an integer, then it is likely that the T-function reflection 
formula T(x)r(l — a;) = vr/ sin(7ra;) will introduce phase 
factors in the residues of poles within C'(—M). This can 
be seen in (110(1 and 112(1 which respectively give phases 
factors of the type (— x) m sin(a7r) and (— x) a+m sin(a7r). 
This leads to terms in the series that will, in general, be 
complex if a is not an integer. Care should be exercised, 
for example, in situations where large-order perturbation 
results are used - since expressing T(3n + 1/2), say, as 
r(3rt) for large n could lead to different imaginary parts. 

It is not always the case that applying (JSJ to a non- 
alternating series will give a finite imaginary part - for 
example, consider the geometric series, ^x 11 , for x = 2, 
say, which when summed gives 1/(1 — 2) = —1. 

An interesting feature of 1(23(1 is that the divergent se- 
ries representation ((25(1 . after the application of © yields 
two convergent series, 1(29(1 and 1(30(1 . which are also valid 
when 9 < 0, even though the integral in ((23(1 is unde- 
fined. Thus, in some cases, the range of applicability of 
the divergent (and convergent) series exceeds that of the 
integral representation of the function in l(23() . 

In general, the difficult part of a problem is to find the 
exact coefficient functions, and often one has to resort to 
some approximation. The results obtained from © are 
only as good as our ability to find the exact coefficient 
function. The form of © is well suited for analytical 
and numerical evaluation of series via contour integration 
(using packages like Mathematica). By expressing some 
divergent series as a convergent series, or even vice- versa, 



one might be able gain a further understanding of the 
nature of the solution to a problem. 
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